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GENERAL SOLUTION OF THE COMPLEX 4-EIKONAL EQUATION 
AND THE 'ALGEBRODYNAMICAL' FIELD THEORY 



m 
O 
o 

(N 

> 
O 

in 



Vladimir V. Kassandrov 1 

Department of General Physics, People's Friendship University of Russia, Or djonikidze 3, 117419 Moscow, Russia 

We explicitly demonstrate the existence of twistor and ambitwistor structure for 4-dimensional complex eikonal 
equation (CEE) and present its general solution consisting of two different classes. For both, every solution can 
be obtained from a generating twistor function in a purely algebraic way, via the procedure similar to that used in 
Kerr theorem for shear-free null congruences. Bounded singularities of eikonal or its gradient define some particle- 
like objects with nontrivial characteristics and dynamics. Example of a new static solution to CEE with a ring-like 
singularity is presented, and general principles of algebraic field theory ('algebrodynamics') related to CEE are briefly 
discussed. 
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1. Ambitwistor structure of complex 
4-eikonal equation 

The four-dimensional eikonal equation in flat Minkowski 
space-time M (we choose the units where c = 1 ) 



\vs\ 2 



0. 



(1) 



is fundamental nonlinear Lorentz-invariant equation 
which governs the propagation of wavefronts (discon- 
tinuties of fields) in every relativistic theory Efl G] • Its 
complex generalization, with S(x,y, z,t) G C, natu- 
rally arise in geometrical optics and quasiclassical limit 
of quantum theory 0] as well as in GTR Sin- 
gular solutions of eikonal equation and its relations to 
other fields have been studied by H. Bateman Pj and, 
recently, by E.T. Newman et al. 0. Possible inter- 
pretation of singularities as particles was discussed by 
Yu.P. Pyit'ev (for 5D-real eikonal) and A.M. Vino- 
gradov [H] (from general viewpoint); see also our work 

m 

In the version of biquaternionic analysis proposed 
in our works El EH E21 complex eikonal equation 
(CEE) (TIJ plays the role similar to that of the linear 
Laplace equation in complex analysis. Its nonlinearity 
arises as a direct consequence of the noncommutativity 
of biquaternion algebra. 

The method to obtain the general integral of eikonal 
equation, i.e. the solution dependent on an arbitrary 
function (starting from the so called "complete inte- 
gral"), is well known (see, for example, ^HllSj)- We 
not only extend these results to the complex holomor- 
phic case. Making use of the obtained twistor and am- 
bitwistor structure of CEE, we present its general solu- 
tion, in a manifestly Lorentz and gauge invariant form, 
and demonstrate that arbitrary CEE' solution belongs 
to one of two and only two different classes. The proce- 
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dure is extremely simple and can be regarded as gener- 
alization of the known Kerr theorem for shear-free null 
geodesies [141 IT5] . The problem of classification and 
evolution of singularities of the eikonal wavefronts and 
of the related caustics also becomes quite transparent in 
the scheme developed below. 

In the C-case static solutions to CEE do exist, the 
known example is given by the Kerr twisting congruence 
with a ring-like singularity (leading to stationary Kerr 
metric in GTR). In section 2 we present another static 
solution to CEE with singularity of a similar ring-like 
structure. 

For one of the two fundamental classes of CEE' so- 
lutions deep connections with solutions of gauge and 
metric fields can be revealed. In particular, for Maxwell 
field related to these solutions of CEE electric charge 
is necessarily a whole multiple of some "unit", elemen- 
tary one. These and other facts form the premises of 
the algebraic field theory - algebrodynamics - which we 
briefly discuss in section 3. 

To start with, let us represent CEE iTjQ in equivalent 
Pfaffian form 



dS = tpdXip, 



(2) 



where the two SL(2, C)-spinors ip — {if a} and ip — 
{ipA'} are introduced, and dX is the differential of her- 
mitian matrix of space-time coordinates X = {X AA }, 
A, A',... = 0,1. For the latters we'll use the stan- 
dard representation X oa = u = t + z, X 11 = v = 
t — z, X 01 = w — x — iy, X 10 — w — x + iy with 
{x, y, z, t} being the Cartesian coordinates and the time 
respectively. Condition (j2j in components reads 



dS 



d A A'S = fAlpA' 



Q X AA> 

and is equivalent to CEE (TIJ since we have 
det \\d AA ,S\\ = d u Sd v S - d w St\jS = 0, 



(3) 



(4) 
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because of the factor structure of the r.h.s. of Eq.©. 
Conversely, the complex gradient 4- vector Oaa'S should 
be null in view of CEE and, therefore, can be always 
represented in the form (EJ with some spinors ip, ip . 
The latters are projective in nature being defined up 
to a complex scalar multiplier X(X) , in view of the 
following symmetry of Eq. ((2j) : 



ip i > Xip, ip A tp, S S. 



(5) 



We proceed now with equivalent transformation of 
Eq. Q of the form 

dS = pd(XiP) - (pX)diP = (pdp - -/dip = 

= p A d(3 A - 1 A 'd^ A/ . (6) 

Here the pair of 2-spinors W = {ip,(3{ = {ip,Xip} re- 
lated to the space-time points (precisely, to the lightcone 
of the point X ) by the incidence condition 1 



/? = A> ((3 A = X AA i, A .) 



(7) 



forms the projective null twistor W. Analogously, the 
pair of 2-spinors W = {ip, 7} = {ip, ipX} forms the dual 
twistor for which the spinors are related by 



l = pX ( 7 A ' = ip A X AA ) 



r AA ' "I 



(8) 



Moreover, twistors W, W are orthogonal to each other 
because of the identity 

ip/3 ee ip(Xip) = (ipX)ip ee jtp => 

=> < W, W >= ip(3 — 'yip ee 0, (9) 

where by <, > the (redifined, see the footnote above) 
twistor scalar product is denoted. 

A pair of mutually orthogonal (projective) twistors 
is known as an ambitwistor |16II14| . Ambit wistor struc- 
ture is closely related to the flow of null geodesies on a 
generic manifold and, by this, to its Riemannian metric 
structure P2J. In invariant twistor form CEE can 
be written as 



dS = < W, dW > 



< dW,W > 



(10) 



Thus, we see that every solution S(X) to CEE 
may be considered to depend on the space-time co- 
ordinates X only "implicitly", i.e. via the compo- 
nents of projective twistor W(X) = {ip(X), (3(X)} 
defined by the structure of null 4-gradient ©, i.e. 
S = S(ip(X), P(X)) . Below, it's natural to distinguish 
two different classes of solutions for which three compo- 
nents of projective twistor W = {ip, (3} are functionally 
dependent and not. 



'We have no need to use here the accepted multiplication by 
i = V— 1 in the r.h.s. of Eqs.JTJ. This is allowed under the proper 
redefinition of the standard twistor scalar product, see below 



2. First class of solutions to CEE 

Let us assume that for some solution of CEE three pro- 
jective components of related twistor W 6 CP 3 are 
independent with respect to space-time coordinates X 
and in some region of M . Making use of the symme- 
try we can consider then that all four twistor com- 
ponents {ipo' ,ipx> , P , /3 1 } are functionally independent. 
Then in view of the eikonal function S — S{tp,(3) 
will be a function of twistor coordinates while the com- 
ponents of the dual twistor {ip, 7} need to have the form 



dS 



yA 



A' 



dS 



dp A ' ' dil> A ' ' 

Applying now the incidence relation JSJ we get 



(11) 



OS X AA' = _9S 



and, finally, 



P A \iP,X) 



d(3 



dS 

dip A , 



dip/ 



dS dS d(3 A 



dip 



d(3 A dip, 



(12) 



Two Eas. <ll2H implicitly define two unknown spinor 
components ip A > at almost every space-time point X 
(exept at a subset of zero measure, see below). Sub- 
sequently resolving algebraic system l(T2|l for different 
'parameters' X and selecting a continious branch of 
roots we come to some field distribution for the spinor 
ip{X) and, accordingly, for (3{X) — Xip(X) . Substitut- 
ing them into S(ip, (3) we determine the dependence of 
generating function S{ip{X), [3{X)) on space-time coor- 
dinates. 

To make sure that the function S(ip A > , X AA ip A ') 
which we've obtained really satisfy CEE let differentiate 
it with respect to the coordinates X AA and get 

Oaa>S = d AA 'ipB' + -^b(X bb d AA >ipB> + 

+ 5*% ip B> ) = P B (d AA ,ip B ,) + QpZ^A'- (13) 

In the region of regularity of derivatives d AA npB' 
where they all are bounded, we have in account of 
Eqs.JTJJ: 



OaA'S = -i^-lpA' = <PAlpA> , 

o/3a 



(14) 



so that the eikonal equation det H^aa'^H = is identi- 
cally satisfied. 

On the other hand, the locus where the derivatives of 
ip(X) are singular corresponds to the branching points of 
C -valued solutions (multiple roots) of algebraic system 
ft2]). i.e. to the caustic condition 



D = det I 



d 2 S 



dip a 1 dip 1 



= 0. 



(15) 
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Thus, we have proved the following 

Theorem 1. Every (analytical) solution S(X) to 
CEE for which all three components of related projective 
twistor W = {ip(X), Xtp(X)} are functionally indepen- 
dent, can be generated by some arbitrary (holomorphic) 
twistor function S(ip,Xtp) via its differentiation and 
subsequent resolving of algebraic system with re- 

spect to ipA' ■ Such solutions have branching points 
determined by condition fifty . Conversely, any function 
S(ip, Xip) for which conditions flSfy can be resolved, 
define a solution of CEE. ■ 

Of course, the procedure can be simplified by means 
of the choice of a special gauge. Making use of partial 
arbitrariness of the spinor ip implied by the symmetry 
JU we can reduce to unity one of the components of the 
spinor ip , say Vo' • Then twistor components acquire 
the form 

V» ' = 1, ipv = G, 0° = wG+u, P 1 = vG+w, (16) 
and generating twistor function reduces to 



S = S^ 1 ',^, I3 1 ) = S{G, wG + u,vG + w); 



(17) 



it should be then differentiated with respect to G, and 
algebraic equation 



(18) 



in the region of regularity, results in some solution G(X) 
and, subsequently, - in some solution S(G(X), (3(X)) of 
CEE. Singularity condition 111 -"it reduces by this to 



D 



d 2 S 
dG 2 



0. 



(19) 



In the form presented by I|17I1 . II18|1 . H19JI the procedure 
becomes much similar to that used for determination 
of singularities of caustics and wavefronts in the theory 
of differentiable mappings and to the method 

discribed in recent works of E.T. Newman with his 
colleagues jS]. However, twistor methods make the 
procedure quite transparent and preserve Lorentz and 
gauge invariance of CEE. We'll further see also that for 
S{X) £ C another class of CEE' solutions does exist 
which has no analogues in the real case. 

Even for the class of solutions dealt with above, un- 
like the real case where always the light-like evolution 
of wavefronts always takes place, a new type of static 
solutions (for which both G(X) and S(X) don't depend 
on time at all) can be obtained in the C-case since the 
equation 



!V>'^(g) 2 + (§ 2 + (g 2 = o 



dy 



dz ' 



can then possess nontrivial solutions. As an example, in 
order to obtain a static axisymmetrical solution, let us 
take the generating twistor function 



S 



G 2 



G 2 



Gf3° - (3 1 + 2iaG wG 2 + 2z*G - w 



(20) 



where constant positive a <E R and z* = z + ia. Apply- 
ing the procedure described above we easily find 



G= ™ 



and for CEE' solution correspondent to i120ll 



S = 



ly 



ww + z* x 2 + y 2 + (z + ia) 



2 ' 



(21) 



(22) 



Both functions jSH i jSH are single- valued, and the 
spinor field G is regular on the whole 3-space. How- 
ever, the eikonal function S as well as its complex (null) 
gradient vector VS become singular at a ring of radius 
'a', i.e. at 



z = 0, x + y 



2 -a 2 . 



(23) 



The solution H22H has, therefore, much in common with 
the famous Appel-Kerr static (stationary) solution of 
the equations of shear-free null congruence (see below, 
section 3) whose singular locus is of the same ring-like 
structure as in H23H . However, in the case considered 
we have a great number of static (axisymmetrical) solu- 
tions to CEE for which the singular locus is bounded in 
space; we expect to study them in detail, together with 
nonstatic solutions, in the forthcoming paper. 



3. Second class of CEE' solutions and 
related physical fields 



It remains for us to consider the case when three com- 
ponents of projective twistor arising in the form © of 
CEE are functionally dependent 2 . Making again use 
of the projective structure of spinors ip, <p manifested 
by the symmetry (0 we can choose the spinor ipA' in 
such a way that any three components of the full twistor 
W = {tp, Xip} would be functionally dependent. This 
means that in the case considered, for every solution of 
CEE, two algebraic constrains between twistor compo- 
nents should exist, of the form 



n (c) (w) = n 



= TT(C) 



{ipA' , X AA ' Tp/ 



= o 



(24) 



where U^- c \ C = 0, 1 are two arbitrary and indepen- 
dent (holomorphic) function of four twistor 'coordi- 
nates'. Analogously to Eas. (jl2fl from the previous sec- 
tion, algebraic system Il24t in implicit form can define, 

2 For a generic twistor, two independent components do always 
exist, see e.g. 1231 
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for a pair of n( c ' being choosed and in the region of 
regularity, a solution for the spinor ip(X) . 

Let some solution ipA'(X) of Ens. <l24H and, conse- 
quently, (3 A (X) = X AA il) A '{X) is found. Then m view 
of the form © of CEE every (analytical) twistor func- 
tion S(W(X)) = S(i/>(X),p(X)) can in principle be a 
solution to CEE (though really it would depend only on 
two functionally independent twistor arguments). 

To make it sure we first differentiate with respect to 
coordinates X AA the generating system J24II and get 



9a A' ipi 



' d(3 A ' 



■(C)B' 



^B'A^l 



(25) 



where by <&b'A the quantities in square brackets are 

1 . 



denoted, and Qt C \ B > is the matrix inverse of 



)(C)B' 



dips 1 



(26) 



Note that in view of Eqs.J2HJ solutions ipA'(X) of 
algebraic system 12411 are defined exept in the singular 
locus represented by the condition 



Q = det ||Q 



(C)B'i 



det 



dyj B > 



= 0. 



(27) 



We'll see further that at these points which correspond 
to multiple roots of Eos. ll24|) the strengths of associated 
physical fields become singular. 

Let us now make use of general expression 11 -ill for 
derivatives of a twistor function S(ip(X), fl(X)) . Sub- 
stituting into it the current expression jSSt for spinor 
derivatives we get 



Oaa'S 



~P B '<S>B'A + dS/d(3 A i> A > = <PA^A> (28) 



where P B are the same as in En. ltl2|l and the spinor 
ip a is now defined by the expression in square brackets. 
We thus see again that the complex 4-gradient Oaa 1 S is 
identically null because of its factorized structure mani- 
fested by the r.h.s. of Ens. p8]t . so that CEE is satisfied 
for every generating twistor function. In this way we 
come to the following 

Theorem 2. Every (analytical) solution to CEE for 
which three components of related projective twistor W 
are funcionally dependent, can be generated by some two 
arbitrary and independent twistor functions US G \il), Xtp) 
via resolving the algebraic system [2$ with respect to 
ip{X) , and is represented by arbitrary twistor func- 
tions S(ij>(X), Xip(X)) 'built' on the solutions ip(X) . 
Branching points for S(X) and ip{X) are defined by 
the condition \2T\) . Conversely, for any pair of 
for which the constraints [2$ can be resolved, a class of 
solution of CEE S(ip(X),XiJ>(X)) is defined. ■ 

As in the previous section the procedure can be es- 
sentially simplified by fixing of gauge. Let one of the 



components of the main spinor, say tpo' , be nonzero in 
some region of space-time. Then it can be brought to 
unity and we come to the gauge 11 (ill used before. Now, 
however, three projective twistor components llKit are 
considered to be functionally dependent; the constrains 
j2H reduce to one algebraic equation of the form 



IL(ip v , (3° , (3 1 ) = U(G,wG + u,vG + w) = 



(29) 



which implicitly define an analytical solution G(X) ev- 
erywhere exept in the branching points given by the 
reduced form of Eos ^27(1: 



Q 



da 

dG 







(30) 



and, of course, on the curves where G{X) becomes zero 
or infinite. Having available a solution G(X) of Eo. ll29|l 
we may be sure that any generating function of three 
projective twistor components 'built' on G{X) , i.e. 

S = S(G(X),wG(X) + u, vG(X) + w) 

(precisely, of any two functionally independent of them) 
should satisfy CEE. For some solutions to CEE i/jqi can 
be zero; in this case another gauge can be used to obtain 
them in analogous manner. 

It's worthy of notice that solutions G(X) of gener- 
ating algebraic equation H29H define geometrically the 
shear-free null geodesic congruences |14l ITB] which are 
related, in particular, to Riemannian metrics of Kerr- 
Schild type. Indeed, Ens.p^f . in the gauge we have 
used, reduce to the system 

d u G = k, d w G = kG, c\tG = 1, d v G = lG, (31) 

where k(X),l(X) are two remaining nonzero compo- 
nents of the matrix $b'A defined by Eas, l|25JI . Elim- 
inating them from Ecis J31Jl we come to the following 
nonlinear system: 



d w G — Gd u G , d v G — Gc\jjG , 



(32) 



which is known to define the principal spinor of a shear- 
free null congruence and whose general solution, in ac- 
cord with the Kerr theorem , is implicitly given 
by Eq.JUJ. 

As to the gauge-free representation, Eas. 1251 can be 
rewritten in the following invariant matrix form: 



(33) 



($ = {&b'a(X)}) which has been earlier proposed in 
our works OElElEBj) m the framework of biquater- 
nionic generalization of Cauchy- Riemann equations and 
which was interpreted as a basic dynamical system (non- 
linear, over-determined, nonLagrangian) for some pe- 
culiar algebraic field theory ( 'algebrodynamics ' ) . The 
latter deals with nontrivially interacting 2-spinor field 
tp{X) and 2 x 2-matrix (complex 4-vector) field 
and was studied also in [201 1191 12~2"] . 
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Let us briefly enumerate the most important features 
of this theory (for details see [23 E])- For system 1(331) 
an unusual gauge symmetry does exist, the so called 
'weak one', which corresponds to the transformations of 
the form 



1p I * A(W)V>, <5>B'A >-» ®B>A + d A A*X 



(34) 



with a gauge parameter A(W) = X(tp,(3) which is al- 
lowed to depend on coordinates X AA only implicitly, 
via the components of the trasforming spinor *p(X) and 
of its twistor counterpart /3(X) = Xtp(X) . 

Owing to the gradient-wise type of its transforma- 
tions, the matrix field &b'a(X) can be interpreted as 
4-potentials of some (C-valued) 'weak' gauge field. In- 
deed, the integrability conditions of Eos. ll3"3l) ddtp = 
result in the self-duality of (matrix) field strengths cor- 
respondent to potentials &(X) and, in its turn, - in 
fulfilment of gauge fields' equations. Precisely, both C- 
valued Maxwell and Yang-Mills equations are satisfied 
identically on the solutions of Eqs.lJHBJ, for the trace 
and the trace-free part of the full matrix field strength 
respectively. For the spinor of complex electromagnetic 
field strength F^b)(X) a beatiful twistor-like expres- 
sion was obtained in 23 1: 



Ft 



(AB) 



1 

Q 



u A1 



_d_ /UaUb 
dG\ Q 



(35) 



where 11,4,11,43, A, B = 0,1 denote the (1-st and 2-d 
order) derivatives of some function H(G, (3°, /3 1 ) with 
respect to its twistor arguments (3 A , and Q = dH/dG 
is the same as in Eci. ll3()|) , If for any II (after differ- 
entiation) the solution G(X) of generating En. (l29]) is 
substituted into I|35J1 the spinor F(ab) will identically 
satisfy vacuum Maxwell equations (exept in the points 
of singularities given by Q = ). 

Thus, for any solution G(X) of Eci. l|29ll condition 
l|30JI defines the shape and the time evolution of singu- 
larities of some associated electromagnetic field. Such 
singular loci can be zero-, one- or two dimensional. In 
the case when the singularities are bounded in 3-space 
they define some singular, particle-like objects with non- 
trivial dynamics simulating interaction of particles and 
even their transmutations. 

Moreover, the value of electric charge associated with 
bounded singularities of electromagnetic field is nec- 
essarily quantized (i.e., equal to the whole multiple of 
some mimimal charge) PI 1101 l2"l] . Such "elementary" 
charge is carried, in particular, by the source of the ex- 
eptional static, axisymmetrical Appel-Kerr solution 



G 



w 



x + iy 



(36) 



z* ± r ~~ ( z + ia) ± ^Jx 2 +y 2 + (z + iaf 
generated by the function Q29JI of the form 

n = G/3°-/3 1 =wG 2 + 2z*G-w = {). (37) 



Its ring-like singularity (|23l) (branching curve of the 
G(X) ) is the sameas for the first-class solution of CEE 
(12211 : it might be interpreted then as a source of elec- 
tromagnetic field defined by (the real part of) Ecis, l|35j) . 
The electromagnetic field itself is identical to that of the 
field of the Kerr-Newman solution in GTR but can carry 
only 'unit', elementary charge (q = ±1/4 in the ac- 
cepted system of units). General theorem on the charge 
quantization has been proved in [231 - This fact, inci- 
dentally, revives the well-known Carter-Lopez model of 
electron [2^1 ESI ED] In the limit a -> the field turns 
into the 'unit-charged' Coulomb field. 

In conclusion, we remind that Ens. ll33|l which we dis- 
cussed above in the framework of general algebraic field 
theory are actually equivalent to the equations of shear- 
free null geodesic congruences i.V2\) and represent the 
second class of solutions to CEE we studied above. 

4. General solution of complex eikonal 
equation 

We are going now to bring together the results obtained 
in theorems 1 and 2 and to describe the general solution 
to CEE (0|. 

Theorem 3. Let in some open subspace O C M there 
exists an analytical function S{X) £ C which is there 
a solution to CEE. Then the compoments of projec- 
tive spinor ipA'(X) and related projective null twistor 
W(X) = {ijj(X),Xi)(X)} are defined via the struc- 
ture of null 4-gradient Oaa'S, and for every solution 
S = 5(W ) is a twistor function, i.e. depends on coor- 
dinates X only implicitly, as follows: 



S(X) = S(tP(X),XiP(X)). 



(38) 



In l(38jl the spinor function tp(X) satisfy either the sys- 
tem {I2l or the system l|24*j) . 

Conversely, any analytical solution to CEE belongs 
to one of two different classes and can be generated by 
one S'(W) or two n( c ), C = 0, 1 arbitrary and indepen- 
dent twistor functions for which algebraic equations fiijj) 
or {24V can b e (analytically) resolved in some subspace 
O C M with respect to ip(X) . 

In the first case (when dS/dip A — 0), all (projective) 
twistor components are independent, and the solution of 
CEE is represented by the function S(W(X)) itself, the 
branching points being defined by En. lTl5ll . 

In the second case (when U.^ c > =0), only two pro- 
jective twistor components remain functionally indepen- 
dent. Then the solution is represented by arbitrary (an- 
alytical) function S(W(X)) and implicitly depends on 
X through two independent twistor arguments only. 
Branching points are then defined by Eo. ll2?jl . ■ 

To conclude, we show that the Cauchy problem for 
CEE with arbitrary analytical initial data S(x, y, z, t ) = 
Sc being given at the moment t — to can be resolved 
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in a simple algebraic way using the procedure above- 
presented. Indeed, one immegiately obtains at t = to 
the 4-gradient Oaa'Sq = PAipA' and the three com- 
ponents of one of the projective twistors, say W = 
{ip,Xip} (in the matrix X the time coordinate is also 
settled to to)- Now it's a trivial step to find out if 
these three components are functionally independent or 
not. In the second case, we immegiately obtain (in the 
gauge for the spinor used) the constrain condition (12911 
represented by a twistor function II(W) . Similarly, in 
the first case we algebraically express the spatial coordi- 
nates x, y, z via three independent twistor components 
and after substitution into Sc come again to a twistor 
function l|17l) >5(W) . Then, in both cases, disposing 
of generating twistor functions II or S and consider- 
ing now the time coordinate in X to be free we are 
able to obtain the full time-dependent solution to CEE 
corresponding to the given Cauchy data by the prose- 
dure described above (i.e. by solving respectively the 
equations II = or dS/dG = with respect to the pro- 
jective spinor component G{X)). The examples of this 
quite transparent procedure will be presented elsewhere. 
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